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Abstract:  The  study  examines  the  capability  for  reduction  of  the  resonant  pressure  fluctuations  in  supersonic  flow 
over  an  open  rectangular  cavity  using  repetitively  pulsed  energy  deposition,  computationally.  The  simulations  are 
performed  for  the  cases  both  with  and  without  energy  deposition  for  the  unsteady  flow  over  a  three-dimensional 
cavity  at  a  free  stream  Mach  number  of  1.5.  The  Reynolds-averaged  compressible  time  dependent  Navier  Stokes 
equations  in  three  dimensions  are  solved  using  the  flow  solver  GASPex.  The  k-<o  turbulence  model  is  used.  A 
Gaussian  temperature  profile  is  utilized  to  model  the  energy  pulse  assuming  the  density  within  the  spot  is  initially 
uniform.  The  results  indicate  that  energy  deposition  reduces  the  pressure  fluctuations  changing  the  flow  structure 
within  the  cavity. 
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NOMENCLATURE 


D:  Cavity  depth 

E:  Total  energy  deposited 

L:  Length  of  cavity 

M:  Mach  number 

Re;  Reynolds  number 

SPL:  Sound  pressure  level 

U:  Free-stream  velocity 

W;  Width  of  the  cavity 

At:  Time  step 

fm:  Rossiter  frequency 

tR  Rossiter  period 

tf:  Final  time 

t,:  Initial  time 

8:  Boundary  layer  thickness 

e:  Non-dimensional  energy  deposition 

parameter 


INTRODUCTION 

High  speed  flows  over  open  cavities  can  produce 
complex  unsteady  tlowfields  that  are  physically 
important  and  become  a  practical  concern  in  aerospace 
applications.  These  complex  unsteady  tlowfields  include 
not  only  the  small-scale  pressure  fluctuations  typical  of 
turbulent  shear  flows  but  also  a  significant  resonance, 
the  frequency  and  amplitude  of  which  depend  on  the 
cavity  geometrj’  and  external  flow  properties. 

Supersonic  flow  past  a  cavity  has  numerous  applications 
in  store  carriage  and  release.  Internal  carriage  of  stores  is 
used  for  supersonic  aircraft  in  order  to  minimize 
aerodynamic  drag,  aerodynamic  heating  and  radar 
signature.  It  also  provides  several  advantages  including 
enhanced  maneuverability  and  expanded  flight  envelope. 
However,  the  intense  pressure  fluctuations  and 
significant  resonant  acoustic  modes  that  are  generated  by 
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the  flow  past  an  open  cavity  at  supersonic  speeds  can 
damage  the  structure  of  the  aircraft  and  stores  and 
impede  successful  store  release.  This  severe  acoustic 
environment  within  the  cavity  can  also  represent  a 
potential  hazard  to  sensitive  instrumentation. 

The  flowfield  structure  within  a  rectangular  cavity  can 
be  categorized  into  two  main  types,  namely  open  and 
closed  cavity  configurations  depending  on  the  cavity 
streamwise  length  to  depth  ratio  LTD  as  illustrated  in  Fig. 
1.  The  flow  in  the  shallow  cavity  where  the  depth  to 
length  ratio  is  small  is  called  closed  cavity  flow,  whereas 
the  flow  for  the  deep  cavity  is  called  open  cavity  flow. 
In  the  open  cavity,  the  free  shear  layer  attaches  on  the 
rear  face,  thereby  forming  a  single  recirculation  region. 
In  the  closed  cavity,  the  shear  layer  reattaches  on  the 
cavity  floor,  and  the  flowfield  is  a  combination  of  a 
backward  and  forward  facing  step. 
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Figure  I.  Types  of  Cavity  Flows 

Numerous  research  studies  on  supersonic  cavity  flow 
have  been  conducted.  A  review  of  the  early  work  in 
subsonic  and  supersonic  flow  past  a  cavity  was  made  by 
Charwat  et  al  (1961).  Rockwell  and  Naudaschcr  (1978) 
reviewed  the  field  and  categorized  flow  types.  Pemg 
(1996)  presented  a  more  recent  review.  Rossiter  (1964) 
proposed  an  empirical  formula  for  predicting  the 
frequencies  of  the  pressure  oscillation  modes,  based  on 
his  shadowgraph  studies  of  subsonic  and  transonic 
shallow  cavities.  Modes  can  be  predicted  well  by 
improved  versions  of  his  formula,  but  any  simple 
method  does  not  exist  for  predicting  the  fluctuation 
amplitudes.  Heller  and  Bliss  (1975)  modified  Rossiter's 
formula  for  the  prediction  of  mode  frequencies.  Bauer 
and  Dix  (1991)  present  a  simple  analytical  method, 
which  uses  the  modified  Rossiter  equation  to  predict 
frequencies. 

Several  experimental  and  computational  studies  have 
also  been  performed  on  cavity  flow.  Leu  and  Dolling 


(1997)  conducted  a  series  of  tests  on  high  speed, 
rectangular,  open  cavity  flows  with  store  release  and 
shock  impingment.  They  measured  fluctuating  surface 
pressure  in  a  cavity,  where  L/H=3  at  Mach  5.  Their 
results  show  that  the  cavity  flow  field  remains 
essentially  the  same  whether  a  store  is  placed  inside  the 
cavity  or  not.  Unalmis  et  al  (1999,  1998)  made  planar 
laser  imaging  and  fluctuating  pressure  measurements  in 
a  Mach  5  cavity  flow,  with  varying  L/D.  They  concluded 
that  the  acoustical  phenomena  inside  the  cavity  do  not 
have  much  effect  on  the  shear  layer  dynamics,  which 
suggests  a  weak  shear  layer/acoustics  coupling  in  the 
current  cavity  flow,  resulting  in  purely  acoustic  modes. 
Hankey  and  Shang  (1979)  predicted  the  pressure 
oscillations  for  supersonic  flow  over  an  open  cavity  by 
numerically  solving  the  two  dimensional  unsteady 
Navier  Stokes  equations  for  Mach  1.5  configuration  of 
Heller  and  Bliss  (1975)  and  obtained  good  agreement  for 
both  the  magnitude  and  frequency  of  the  lowest  order 
modes.  Rizzetta  (1988)  performed  a  three  dimensional 
unsteady  RANS  simulation  at  a  free  stream  Mach 
number  of  1.5  and  observed  good  agreement  with 
experimental  data  with  which  comparisons  were  made  in 
terms  of  mean  static  pressure  and  overall  acoustic  sound 
pressure  levels  within  the  cavity.  Shih  et  al  (1994) 
developed  a  numerical  procedure  for  the  simultaneous 
implicit  solution  of  the  coupled  k-£  and  Navier-Stokes 
equations  for  compressible  viscous  flows  and  compared 
the  results  with  experiments  and  other  computations. 

A  variety  of  active  and  passive  control  techniques  have 
been  investigated  to  reduce  the  sound  pressure  level 
(SPL)  associated  with  the  resonant  acoustic  fluctuations. 
In  passive  control  approach,  the  cavity  is  modified  in 
some  way,  such  as  a  mounted  device.  Passive  control 
methods  are  inexpensive,  simple  and  effective  in 
suppressing  the  oscillations.  However,  the  performance 
of  the  cavity  in  time  varying  conditions  may  be  worse 
than  that  without  control,  since  passive  control  uses 
permanent  devices.  Active  control  methods  can 
continuously  change  to  adapt  to  varying  flow  conditions. 
Examples  of  active  and  passive  control  are  the  studies 
performed  by  Lamp  and  Chokani  (1997)  and  Stallings  et 
al  (1994).  Lamp  and  Chokani  (1997)  performed 
computations  on  cavity  flows  with  suppression  using  jet 
blowing.  Stallings  et  al  (1994)  investigated  the  effect  of 
passive  venting  on  static  pressure  distributions  in 
cavities.  Zheltovodov  et  al  (1991)  also  studied 
supersonic  flow  over  rectangular  cavities  and  methods  of 
flow  control.  Aradag  et  al  (2004)  studied  the  effects  of 
laser  energy  deposition  on  supersonic  cavities 
computationally  and  Lazar  et  al  (2008)  studied  the  same 
cavity  configuration  experimentally. 

OBJECTIVES 

The  aim  of  this  research  is  to  simulate  supersonic 
turbulent  flow  past  an  open  rectangular  cavity,  and  to 
examine  the  capability  of  pulsed  energy  deposition  to 
reduce  the  resonant  pressure  levels  in  the  cavity. 
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STATEMENT  OF  THE  PROBLEM 

The  coordinate  system  and  the  configuration  that  are 
considered  are  given  in  Fig  2.  The  geometry  represents 
one  of  the  cases  of  the  experimental  data  of  Kaulman  et 
al  (1983).  The  length  to  depth  ratio  for  the  cavity,  L/D  is 
5.07,  the  length  to  width  ratio,  L/W  is  1.90.  The  other 
principal  parameters  are  the  free  stream  Mach  number  of 
1.5  and  the  Reynolds  number  of  1.09  10*  based  on  the 
free  stream  conditions  and  the  cavity  length,  L.  Tlie  free 
stream  static  temperature  is  218  K  and  the  free  stream 
total  pressure  is  66.4  kPa.  The  ratio  of  the  boundary 
layer  thickness  immediately  upstream  of  the  separation 
point  to  the  cavity  depth  5/D  is  0.213. 


A  periodic  pulsed  energy  input  of  1 .0  mJ  is  given  to  the 
leading  edge  of  the  cavity  to  simulate  the  energy 
deposition  by  laser,  which  corresponds  to  a 
dimensionless  energy  parameter,  e=4.0.  non- 
dimensionalized  using  static  enthalphy.  The 
dimensionless  energy'  deposition  parameter  (Knight  et  al. 
2003)  is  defined  as 


The  numerator  is  the  amount  of  energy  deposition  in  V 
in  time  Te.  The  denominator  is  the  sum  of  two  terms. 
The  first  term  is  the  static  enthalphy  in  region  V  at  the 
beginning  of  the  energy  pulse.  The  second  term  is  the 
static  enthalphy  flux  through  the  cross  sectional  area  A 
during  the  characteristic  timescale  of  the  energy 
deposition  Tj  (Knight  et  al,  2003).  The  second  term  is 
neglected  in  this  study  since  the  energy  pulse  is  assumed 
instantaneous. 

A  spherically  symmetric  initial  temperature  distribution 
is  proposed  to  model  each  energy  pulse  (e.g.,  by  laser) 
assuming  the  energy  is  added  instantaneously  at  constant 
volume  (Yan  et  al,  2003).  Therefore  the  density  is 
constant  and  an  ideal  gas  is  assumed.  The  temperature 
variation  using  a  Gaussian  profile  can  be  written  as: 

AT  = 

(2) 


where  the  peak  temperature  variation,  ATo  is  determined 
by  the  total  energy  de|X)sited,  E,  which  is  1.0  mJ.  The 
energy  deposited  is  given  by 


E  = 


.sin 


(3) 


where  Cv  is  the  specific  heat  at  constant  volume. 
Substituting  (2)  into  (3)  and  integrating,  we  obtain 


Aro  = 


(4) 


where  To  is  related  to  the  assumed  initial  radius  Ro=0.9 
mm  of  the  energy  pulse  obtained  from  the  perturbation 
focal  volume  Vo  =  4/37tRo^  and  set  to  be  Ro/2.  From  Eq. 
(2),  AT  will  reach  2%  of  ATo  at  r=Ro.  The  ratio  of  the 
energy  pulse  radius  to  incoming  boundary  layer 
thickness  Ro/6  is  0. 1 77. 

Despite  its  inherent  simplicity,  the  energy  deposition 
model  provides  an  accurate  prediction  of  the  flowfield 
generated  by  a  laser  pulse  outside  of  the  laser  focal 
volume  and  its  immediate  neighborhood.  The  model  was 
validated  by  comparison  with  Filtered  Rayleigh 
Scattering  measurements  (Yan  et  al.  2003). 

NUMERICAL  METHODOLOGY 


The  flow  solver  is  GASPex  Version  4.1.0+  (1996)  which 
is  a  structured  multiblock  CFD  solver  that  is  applicable 
to  compressible  flow  fields  approximately  at  Mach  0.2 
or  higher.  It  solves  the  Reynolds-averaged  compressible 
time  dependent  Navier  Stokes  equations  in  three 
dimensions,  utilizing  a  finite  volume  spatial 
discretization  in  which  the  state  variables  are  stored  at 
the  cell  centers.  Inviscid  fluxes  are  modeled  using  third 
order  upwind  Van  Leer  scheme.  The  computations  are 
second  order  accurate  both  in  time  and  space. 

The  turbulence  model  is  k-O),  which  is  a  two-equation 
model.  This  turbulence  model  is  widely  used  to  simulate 
incompressible  and  compressible  flows.  In  this  model, 
transport  equations  for  the  turbulent  kinetic  energy  k  and 
a  second  parameter  to  (the  rate  of  dissipation  of  energy 
per  unit  volume  and  time)  are  solved  (Wilcox,  1993). 
The  Sutherland  viscosity  law  is  employed  together  with 
the  ideal  gas  relation  in  the  computations. 

The  configuration  is  symmetric  about  the  cavity 
centerline,  so  the  computational  domain  is  half  of  the 
physical  domain.  Inflow  boundary  conditions  were 
obtained  from  the  numerical  solution  of  the  two- 
dimensional  steady  flow  equations  over  a  flat  plate  at  the 
free  stream  conditions  of  the  cavity  flow,  using  program 
EDDYBL  (Wilcox,  1993).  This  computation  employed  a 
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finer  grid  distribution  than  that  is  used  for  the  cavity 
computation,  so  a  cubic  spline  interpolation  was  used  to 
interpolate  the  results  to  the  cavity  flow  computational 
domain.  The  flat  plate  computational  domain  extended 
from  the  leading  edge  of  the  cavity  experimental 
configuration  (x  =  -1.526L)  to  the  forward  bulkhead. 
Results  of  this  solution  at  x  =  -0.249L  established 
upstream  profiles  of  the  dependent  variables  at  all 
spanwise  locations.  The  incoming  boundary  layer 
thickness  matched  the  experiment.  On  the  solid 
boundaries,  the  no-slip  conditions 


from  the  cavity  symmetry  plane.  Stretching  takes  place 
along  both  sides  starting  from  the  centers  of  the  y-z 
planes  of  the  computational  domain  up  to  27  z-grid 
points.  ICEMCFD  software  was  used  for  the  generation 
of  structured  grid.  At  the  inflow  boundary,  there  are 
approximately  24  points  within  the  boundary  layer.  The 
extent  of  the  simulation  domain  was  -0.24865<x/L<  1 .5, 
-1 .0<  y/D<4.0, 0<z/(W/2)<2.0. 

The  unsteady  solutions  were  obtained  using  the  implicit 
dual  step  method.  The  time  increment.  At  is  Ips 
corresponding  to  0.0012  non-dimensional  time,  non- 
dimensionalized  by  the  fundamental  Rossiter  period,  (tn) 
The  Rossiter  period  is  the  inverse  of  Rossiter  frequency 
which  is  defined  as 


was  employed  along  with 


(6) 


T  =  r„„,, 


(7) 


where  the  isothermal  wall  temperature  was  taken  as  the 
adiabatic  wall  temperature  that  was  calculated  to  be 
304.8  K  and  n  is  the  direction  normal  to  the  surface.  The 
solid  wall  boundary  conditions  for  k  and  to  arc 


k  =  0  (8) 

fnilr 
u;  = - 

t^l  (9) 

where  U*  is  the  wall  shear  velocity,  p  is  the  local 
density,  p  is  the  local  laminar  viscosity  and  r  is  a 
parameter  related  to  the  non-dimensional  roughness.  For 
the  outflow  and  upper  computational  boundaries,  first 
order  extrapolation  was  used.  On  the  centerline, 
symmetry  conditions  were  satisfied.  Above  the  cavity 
mouth,  the  initial  conditions  were  taken  as  the  inflow 
conditions.  Inside  the  cavity,  wall  values  were  assigned. 


/«  = 


mU 


(10) 


where  p  is  the  number  of  waves  or  mode  number,  L  is 
the  cavity  length,  LI  is  the  free  stream  velocity.  Mo,  the 
ratio  of  free  stream  velocity  to  the  stagnation  speed  of 
sound  which  is  given  by  ao  =/  R  To,  is  calculated  to  be 
1.246.  K  is  taken  as  0.55  for  this  cavity  configuration 
(Hankey  and  Shang.  1979). 

For  the  computations,  one  of  the  Beowulf  computer 
clusters  of  Rutgers  University  (mphase)  was  used.  The 
core  of  the  mphase  cluster  consists  of  60  dual  and  single 
processor  Linux  machines.  A  total  of  32  processors  in 
this  cluster  were  used  for  the  computations.  Each  time 
step  takes  148  seconds  with  32  processors.  Each 
processor  uses  500  MB  of  memory.  The  CPU  time  per 
time  step  per  cell  is  6miliseconds. 


RESULTS 


The  simulations  for  unsteady  supersonic  flow  past  an 
open  rectangular  cavity  with  and  without  energy 
deposition  were  conducted  for  one  of  the  Wst  cases  of 
the  experimental  study  by  Kaufman  et  al.  and  the  results 
were  compared  with  the  experimental  data  along  the 
cavity  centerplanc,  for  the  case  without  energy 
deposition. 


The  computational  domain  has  a  nonuniform  Cartesian 
mesh  consisted  of  two  zones  that  are  the  inside  and 
outside  of  the  cavity.  The  total  number  of  cells  in  the 
inside  zone  are  95  X  54  X  35,  and  1 77  X  62  X  70  in  the 
outside  zone,  in  the  x,  y.  z  directions,  respectively.  Grid 
clustering  was  employed  at  several  places  in  the  domain. 
Exponential  stretching  was  used  for  this  purpose  with 
Anmin^  of  1.0,  at  the  solid  boundaries.  The  minimum 
grid  spacing  employed  is:  AXmm 

/L=Aym,„/L=Az,„„/L=0.00004.  Both  within  and  outside 
the  cavity,  the  y  value  is  stretched  from  the  minimum 
value  at  the  mouth  of  the  cavity  to  a  constant  value  of 
0.02.  over  the  first  27  grid  points.  The  first  12  z-grid 
points  have  a  uniform  spacing  of  Az/L=0.02,  starting 


The  flowfield  solution  was  integrated  in  time  starting 
from  the  initial  profiles.  After  1000  ps,  it  was  judged 
that  the  flow  field  was  free  of  initial  transients  and  had 
reached  a  periodic  state  at  which  the  oscillations  were 
self-sustaining.  This  time  corresponds  to  1.2 
fundamental  Rossiter  periods.  Then,  computations  were 
performed  for  an  additional  12  Rossiter  periods, 
corresponding  to  10,000  time  steps,  to  compare  the 
computed  fundamental  Rossiter  period  with  the 
theoretical  value. 

The  pressure  fluctuations  at  a  location,  on  the 
downstream  face  at  the  centerplane.  at  y/D=0.6  where  y 
is  measured  from  the  cavity  floor,  arc  shown  in  Fig  3 
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Even  though  the  mean  pressure  value  is  not  much  higher 
than  the  free  stream  static  pressure,  peak  pressures  can 
be  even  2.5-3  times  the  freestream  value. 


Figure  3.  Pressure  Fluctuations  on  centcrplanc  at  y/D=-0.4 

In  Fig.  4,  the  acoustic  streamwise  pressure  distribution 
on  the  cavity  floor  is  given  as  the  Sound  Pressure  Level 
(SPL).  in  decibels,  defined  as; 
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and  q  is  the  sound  pressure  reference  level  of  2  10  ’  Pa 
and  tf-ti  (Final  time-initial  time)=12tR.  The 
computations  are  within  10  dB  of  experimental  results. 

The  first  energy  pulse  of  1 .0  mJ,  which  corresponds  to  a 
dimensionless  energy  parameter  of  4.0,  was  deposited  to 
the  flow  aAer  12  Rossiter  periods.  An  energy  input  of 
same  amount  was  added  at  the  beginning  of  each 
fundamental  Rossiter  period.  The  effects  of  the  amount, 
location  and  frequency  of  energy  deposition  was  not 
tested.  Energy  is  only  deposited  periodically  at  the 
beginning  of  each  Rossiter  period.  Computations  were 
performed  for  12  Rossiter  time  periods  after  the  first 
energy  deposition. 

Pressure  fluctuations  at  the  experimental  measurement 
location  at  the  aft  bulkhead  after  energy  deposition  for 
12  Rossiter  periods  is  shown  in  Fig.  5.  Energy 
deposition  changes  the  pressure  distribution  as  seen  in 
the  figure.  It  decreases  the  peak  pressures  compared  to 
the  no  energy  deposition  case  (Fig.  3),  whereas  the  mean 
pressure  remains  almost  the  same. 


Figure  5.  Pressure  Fluctuations  on  centerplane  at  y/D=-0.4 
after  energy  deposition 


Figure  4.  SPL  Distribution  on  Cavity  Floor  (Without  energy 
deposition) 
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Figure  6.  Power  Spectrum 
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The  power  spectrum  versus  frequency  graph  which 
shows  the  fundamental  pressure  frequencies  in  the  flow, 
obtained  from  pressure  histories  is  given  for  the  energy 
deposition  case,  after  12  Rossiter  periods,  on  the  s^e 
figure  with  the  case  before  energy  deposition  in  Fig.  6. 

(13)  The  peak  powers  decrease  after  12  Rossiter  periods  with 
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pulsed  encrg)’  deposition  at  the  beginning  of  each 
Rossiter  period.  The  peak  values  of  the  power  spectrum 
are  at  almost  the  same  frequencies  as  the  case  without 
energy  deposition. 

The  acoustic  streamwise  pressure  distribution  on  the 
cavity  floor  at  the  centerplane,  is  given  as  SPL  in 
decibels  in  Fig.  7  after  energy  deposition  for  1 2  Rossiter 
periods.  The  pressures  decrease  uniformly  along  the 
cavity  floor  after  energy  deposition. 


Figure  7.  Pressure  Disuibution  on  Cavity  Floor  before  and 
after  Energy  Deposition 


In  Fig.  8  and  Fig.  9,  sound  pressure  level  distributions 
are  shown  on  the  forward  and  aft  bulkheads,  at  the 
centerplane.  after  12  Rossiter  periods.  They  also 
decrease  uniformly  after  energy  deposition.  In  Fig.  8.  at 
the  forward  bulkhead,  the  location  where  energy  is 
deposited  remains  the  only  position  where  the  SPL 
values  do  not  decrease,  which  is  reasonable  since  the 
pressures  are  increased  at  the  energy  deposition  location. 


*  y/D 

Figure  8.  Pressure  Distribution  on  Upstream  Face  of  the 
Cavity  before  and  after  Energy  Deposition 

Energy  deposition  is  amtinued  for  an  additional  10 
Rossiter  periods  to  examine  the  convergence  of  the 
power  spectrum  at  the  experimental  measurement 
location.  The  power  spectrum  docs  not  converge  to  a 
single  state  when  the  computations  are  continued  for 
several  Rossiter  periods.  The  peak  powers  of  the  power 


spectrum  for  the  first  and  second  Rossiter  modes  t/tR 
since  the  beginning  of  energy  deposition  are  shown  in 
Fig.  10.  The  peak  power  for  the  fundamental  Rossiter 
frequency  in  the  absence  of  energy  deposition  is  also 
shown.  It  is  clearly  seen  from  Fig.  10  that  the  peaks  of 
the  power  spectrum  change  after  each  Rossiter  period,  if 
we  continue  the  computations  by  depositing  energy  at 
the  beginning  of  each  Rossiter  period.  Although  the  peak 
powers  are  changing  after  deposition  of  energy,  after 
each  Rossiter  period,  they  still  remain  smaller  than  the 
one  before  energy  deposition. 


Figure  9.  Pressure  Distribution  on  Downstream  Face  of  the 
Cavity  before  and  after  Energy  Deposition 


Figure  10.  Peak  Powers  of  the  Power  Spectrum 


CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FUTURE  WORK 

The  capability  for  reducing  the  resonant  pressure 
fluctuations  in  supersonic  flow  over  an  open  cavity  with 
laser  energy  deposition  is  examined.  This  is  the  first 
study  that  gives  us  an  opportunity  to  see  the  effects  of 
energy  deposition  for  the  control  of  pressure  oscillations 
in  the  cavity.  Results  indicate  that  energy  deposition 
reduces  the  pressure  fluctuations  within  the  cavity, 
changing  the  flow  structure  for  several  Rossiter  periods. 

The  phase,  location  and  magnitude  of  energy  deposition 
must  be  optimized  to  have  the  best  results  in  terms  of  the 
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reduction  of  pressure  oscillations.  The  next  step  in  this 
research  must  be  to  see  the  effects  of  phase,  amount, 
location  and  duration  of  energy  deposition  on  the 
structure  of  the  cavity  configuration.  The  limits  for  the 
reduction  of  resonant  pressure  fluctuations  by  energy 
deposition  and  the  way  they  correlate  to  the  parameters 
of  energy  deposition  is  another  concern. 
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